Bicontinuous and mixed gels in binary mixtures of patchy colloidal particles 
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We investigate the thermodynamics and percolation regimes of model binary mixtures of patchy 
colloidal particles. The particles of each species have three sites of two types, one of which promotes 
bonding of particles of the same species while the other promotes bonding of different species. 
We find up to four percolated structures at low temperatures and densities: two gels where only 
one species percolates, a mixed gel where particles of both species percolate but neither species 
percolates separately, and a bicontinuous gel where particles of both species percolate separately 
forming two interconnected networks. The competition between the entropy and the energy of 
bonding drives the stability of the different percolating structures. Appropriate mixtures exhibit 
one or more connectivity transitions between the mixed and bicontinuous gels, as the temperature 
and/or the composition changes. 



I. INTRODUCTION 

Colloidal self-assembly provides a route for the fabrica- 
tion of functional materials with novel mesoscopic struc- 
tures, which are robust but adaptable to changing exter- 
nal conditions. Indeed, soft colloidal solids including gels 
and glasses, are examples of three dimensional structures 
that may be tuned, e.g., by changes in temperature. 

Advances in the synthesis of well-defined colloidal par- 
ticles, with surface patterning in the sub-micron range, 
paves the way for the design of their macroscopic be- 
haviour [1-8]. One focus of research concerns the devel- 
opment of large-scale fabrication techniques required for 
the exploitation of the new materials in a wide range of 
applications [9]. The other, which includes modeling and 
computational efforts, based on the fact that anisotropic 
colloids may be viewed as the molecules of novel materi- 
als, aims at tailoring their (self) assembly into a variety 
of functional structures [l0|, [lil • 

The primitive model of patchy colloids consists of hard- 
spheres with / patches on their surfaces. The patches act 
as bonding sites and promote the formation of well de- 
fined clusters, whose structure and size distribution de- 
pend on the properties of the patches (/ and the bonding 
energy) and on the thermodynamic conditions (density 
and temperature). Computational studies of patchy par- 
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tide models, over the last decade, produced a number 
of results [12|. In particular, Sciortino and co-workers 
showed that, for low values of / (approaching 2), low 
densities and temperatures can be reached without en- 
countering the phase boundary. These empty phases 
were shown to be network (percolated) liquids, suggest- 
ing that, on cooling, patchy particles with low functional- 
ity assemble into amorphous states of arbitrary low den- 
sity - equilibrium gels. These gels are novel structures 
that may be assembled under controlled and reproducible 
equilibrium conditions. 

Remarkably, the results of the simulations are well 
described by classical liquid state theories: Wertheim's 
first order perturbation theory p!3l-[l6| predicts cor- 
rectly the equilibrium thermodynamic properties; Flory- 
Stockmayer p!7l-[T9| theories of polymerization describe 
quantitatively the size distributions of the patchy parti- 
cle clusters as well as the percolation threshold. 

Very recently, Ruzicka and co-workers reported experi- 
mental evidence of empty liquids in dilute suspensions of 
Laponite (lOl confirming the patchy particle model as the 
primitive model of real equilibrium gels. In this paper 
we use theoretical tools to study the structure of equi- 
librium gels in model binary mixtures of patchy colloids. 
We investigate the stability of these amorphous soft-solid 
materials and show that there are four types of gels, the 
stability of which may be controlled by the temperature 
and composition of the mixture. 

The phase diagram of mixtures of patchy particles with 
identical patches was investigated by Bianchi et al [2lj 
and established the stability of the empty fluid regime. 
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Later on we clarified the conditions for the stability of 
empty fluids and characterized their (network) structure 
using a generalization of the Flory-Stockmayer theory 
(22| . The phase diagram of monodisperse systems and 
of mixtures of colloidal particles with dissimilar patches 
was also investigated, elucidating the conditions for the 
emergence of criticality and establishing more general 
conditions for the stability of network (percolated) flu- 
ids dlHiii. 

In the context of a related model, Hall and co-workers 
[261 , [27| investigated the structures formed by monodis- 
perse systems and equimolar mixtures of dipolar colloids 
with particles of two different sizes and dipole moments. 
Their studies, based on (discontinuous) molecular dy- 
namics, focused on kinetic pathways and the resulting 
non-equilibrium structures. In addition to solid phases, 
they report two novel structures, both bicontinuous gels. 
The bicontinuous gels were shown to have tunable pore 
size determined by the size and dipole moment ratios of 
the colloidal particles. 

Here we investigate the thermodynamic phase diagram 
and the percolation regime of a model binary mixture of 
patchy colloidal particles. Each species has three patches 
with two different types of bonding energies, one pro- 
moting bonds between the same species while the other 
promotes bonds between different species. We have iden- 
tified four distinct percolated phases: two gels where only 
one of the species is percolated, a mixed gel where the two 
species are percolated but neither species percolates by 
itself, and a bicontinuous gel where the two species per- 
colate independently, forming two interpenetrating span- 
ning networks. The competition between the entropy of 
bonding and the energy of bonding determines the sta- 
bility of the different network phases. We show that for 
appropriate mixtures a transition between the mixed and 
bicontinuous gels is induced as the temperature and/or 
the composition of the mixture changes. 

The remainder of the paper is organised as follows. 
In section [IT] we describe the model and the theory: 
Wertheim's thermodynamic perturbation theory for mix- 
tures of associating fluids (jll A[) and a generalization of 
the Flory-Stockmayer theory of percolation to mixtures 
of patchy particles (jllBj) . In section [TTIl we present the re- 
sults for a series of representative mixtures. We describe 
the fluid phase diagrams, with emphasis on the percola- 
tion regime, of symmetric and asymmetric mixtures and 
discuss the conditions for the connectivity transitions be- 
tween mixed and bicontinuous gels. Finally, in Sec. IIVI 
we summarize our conclusions. 



II. MODEL AND THEORY 

The model is a binary mixture of Ni and N2 equi- 
sized hard spheres (HSs) with diameter a. Each species 
has three interaction or bonding sites on its surface. 
The association between bonding sites is described us- 
ing Wertheim's first-order perturbation theory. There- 
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FIG. 1. The model: a binary mixtures of equisized hard 
spheres with three bonding sites of two different types on 
each particle. 

fore, two particles can form one single bond between two 
sites, one on each particle. Bonding sites are distributed 
randomly on the particle surfaces in such a way that all 
sites are available to bonding (i.e. there is no shading of 
any site by nearby bonds). 

Species 1 has two bonding sites of type A and one 
of type C while species 2 has two sites of type B and 
one of type D (see Fig. [1]). Although the model may 
seem complex, the idea is quite simple: each species has 
two types of bonding sites, one which promotes bonds 
between particles of the same species (sites of type A 
and B) while the other promotes bonds between different 
species (sites of type C and D). Note that the minimum 
number of bonds required for each species to percolate is 
3. 



A. Helmholtz free energy: Wertheim's 
thermodynamic perturbation theory 

In references p!3l-[l6[ [28| Wertheim's first order pertur- 
bation theory is described in detail. Here we state briefly 
the main results and set the notation. 

The Helmholtz free energy of the binary mixture of 
patchy particles can be written as the sum of two contri- 
butions: 

fH = F/N = fns + /6, (1) 

where fns is the free energy density of the reference sys- 
tem of hard-spheres and the perturbation due to the 
attractive bonding interactions. N = Ni^N2i^ the total 
number of particles. 

The free energy of the reference system of HSs, JhSi 
may be split into ideal-gas and excess terms: Jhs = 
fid + fex ' The ideal-gas free energy is given (exactly) by 

= Inr? - 1 + ^ x(^) ln(x(^)VO, (2) 

i=l,2 

where j3 = kT is the inverse thermal energy, Vi is the 
thermal volume, x^*^ = Ni/N is the molar fraction of 
species i and r] = r]i-\-r]2 '^^ the total packing fraction {r] = 
VsP, with p the total number density and Vs = 7r/6a^ the 
volume of a single particle). The excess part accounts for 
the excluded volume interactions between hard spheres. 
We have approximated it using the Carnahan- Starling 



3 



equation of state for HSs [2^ 
have the same diameter): 



(note that both species 
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(3) 



The bonding free energy, J^, accounts for the interac- 
tions between sites. Within Wertheim's first-order per- 
turbation theory fl3l - [l6l . it is written as [IF 




(4) 

where V{i) is the set of bonding sites or patches on one 
particle of species i {i.e., r(l) = {A, A, C}, r(2) = 
{B^B^D}) and n{T{i)) is the total number of bonding 

sites per particle of species i. The variables {xi'^} are 
the probabilities of finding one site of type a on a particle 
of species i not bonded. The bonding free energy has two 
contributions: the bonding energy and an entropic term 
related to the number of ways of bonding two particles. 

The law of mass action establishes a relation between 
{Xa^} and the thermodynamic variables: 



= 1,2 



(5) 



/S.ai^ characterises the bond between a site a on a par- 
ticle of species i and a site 7 on a particle of species j. 
We model the interaction between sites by square well 
potentials with depths e^^ that depend on the type of 
bonding sites (a and 7) but not on the particle species 
(i and j). As a result, when the particles have the same 

diameter, /S.^^ are independent of the particle species, 
and are given by 



'—^a.'y 
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gHs{r) [exp{/3eaj) - 1] dr. (6) 



9hs{'^) is the radial distribution function of the reference 
HS fiuid and the integral is calculated over the bond vol- 
ume Vb (we have considered that all bonds have the same 
volume, Vb = 0.000332285<j^). We approximate gns by 
its contact contact value (the bond volume is very small) . 
Then Eq. (6) simplifies to: 



= — [exp{f3ea^) - 1] Ao{r]), 



where 



1-^/2 

(1 - 7?)3 



(7) 



(8) 



is the contact value of gHS- By substituting A^^ given 
by Eq. ([7j) in Eq. ([5]) we find that Xa^ depends only on 
a, the type of site (i.e., Xa'' = X^^ V i). 




level i 



level i + 1 




FIG. 2. Schematic representation of a tree tree- like cluster of 
a binary mixture of patchy particles (right). On the left one 
particle of species k (at level i + is bonded to a particle of 
species j (at level i) through sites of type 7 and /3 respectively. 



In what follows we will denote the composition of the 
mixture x by the molar fraction of species 1: x = x^^"^ 
(x^^) = 1 — x). We choose to minimize the Gibbs free 
energy per particle {g = p/ p -\- /h) to obtain the equilib- 
rium properties of the mixture. We set the composition 
X, pressure p and temperature T and locate the binodals 
by a standard common-tangent construction on g{x). 



B. Percolation thresholds 

The percolation threshold is analysed using a gener- 
alization of the Flory-Stockmayer theory of polymeriza- 
tion (see Refs. (22|,[3Q| for details) which neglects 
closed loops. The percolation thresholds and the cluster 
size distributions were found to be in good agreement 
with the results of Monte Carlo simulations for the pure 
fiuids [sH, and we expect the same level of accuracy for 
binary mixtures. 

(k) 

Let denote the number of bonded sites of type 

7 on particles of species k at the level i + 1 of a tree-like 
cluster of particles (see Fig. [2]). This is related to {n--^^}' 
the set of all types of bonded sites in the previous level, 
through the recursive relations given by (see Fig. [2] and 
Refs. [22, 30]): 



{I 



/5 



^(3a ) n, 



(9) 

where: Td{j) is the the set of distinct bonding sites on 
species j {i.e., Td{l) = {A,C} and Td{2) = {B,D}); 
the sum on j runs over all species j = {1,2}; /^-^^ is the 

number of f3 sites on a particle of species j {e.g., f^^^ = 2); 
Sf3a is the Kronecker delta and pf^.^^^ is the probability 
of bonding a site /3 on a particle of species j to a site 7 
on a particle of species k. Summing over all k and 7, we 
obtain 



E E 



(10) 



the probability of finding a bonded site /3 on a particle 
of species j. At this point, it is possible to establish 
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a connection with thermodynamics through the law of 
mass action, since by definition: 

P^,-=l-X«. (11) 

A term- by-term analysis of Eqs. (pT]) and (|5j) yields 
expressions for the equilibrium probabilities pf^.^^^ in 
terms of the thermodynamic variables. 

In order to calculate the percolation threshold it is con- 
venient to write equations ([9]) in matrix form: 

n^=Tno, (12) 

where fii is a vector with components nf}^ and T is a 
square matrix of size VdiX) ^ ^^(2) with: 

T,,a,= Yl (//3'^-^/3c.). (13) 

For the binary mixture of particles with three sites, T 
is a 4 X 4 matrix |32i] : 



Til 


T21 


Tl2 


T22 



^ PAA + PC A ^Paa Pba + PDA ^Pba \ 
^ PAC + Pec '^PAC Pbc + Pdc "^pbc ^-^^^ 
Pab + PCB "^Pab Pbb + Pdb "^Pbb 

\ PAD + PCD '^PAD PBD + PDD "^PBD / 

Let A+ denote the largest (absolute value) eigenvalue of 
T. Then, the progressions defined by Eq. ([T2]) converge 
to if I A+ 1 < 1 and diverge if |A+| > 1. The percola- 
tion threshold occurs when |A+| = 1. The analysis of A+ 
allows us to distinguish percolated from non-percolated 
structures of the mixture, but it does not differentiate, for 
example, bicontinuous from mixed gels as the mixture is 
percolated in both cases. In order to distinguish different 
percolated structures we have to analyse, in addition to 
the percolation threshold of the mixture, the percolation 
thresholds of each species. The matrix T may be split 
into four blocks T^j, i,j = {1,2} (see Eq. (fT4|) ). where 
each block Tij is a 2 x 2 matrix that accounts for (all) the 
bonds between species i and j. The largest eigenvalues of 

~ ~ fll) (22) 

Til and T22 (A^ ^ and A^ ^ respectively) yield the perco- 
lation thresholds of both species. In a percolated binary 
mixture (|A+| > 1) it is then possible to distinguish four 
different percolating structures: 

fii) (22) 

• |Ay < 1| and \X\_ < 1|. The mixture is perco- 
lated but removal of one of the species destroys the 
connectivity of the spanning cluster. We refer to 
this structure as a mixed gel (MG) and illustrate 
it schematically in Fig. [3] (right). 

• |A^^^^ > 1| and |A^^^^ > 1|. The mixture as weh 
as both species are percolated. We refer to this 
structure as a bicontinuous gel or "bigel" (BG) as 



the spanning cluster consists of two interconnected 
spanning clusters of the pure components. Removal 
of one species does not destroy the connectivity of 
the spanning cluster. The bigel is illustrated in the 
middle of Fig. m 

(11) (22) 

• |Ay > 1| and \X\_ ^ < 1|. The mixture and species 

1 are percolated but species 2 is not. Removal of 
species 2 does not destroy the connectivity of the 
spanning cluster. We refer to this structure as a 
standard gel Gi and illustrate it schematically in 
Fig. [3] (left). 

fii) (22) 

• \X\_ < 1| and \X\_ > 1|. The mixture and species 

2 are percolated but species 1 is not. Removal of 
species 1 does not destroy the connectivity of the 
spanning cluster. We refer to this structure as a 
standard gel G2. 

Table |T] classifies the structure of the binary mixture 
based on the analysis of the percolating cluster. 



|A+| 




|Af)| 


percolation state 


< 1 


< 1 


< 1 


non-percolated 


> 1 


< 1 


< 1 


percolated: mixed gel 


> 1 


> 1 


> 1 


percolated: bigel 


> 1 


> 1 


< 1 


percolated: gel 1 


> 1 


< 1 


> 1 


percolated: gel 2 



TABLE 1. Structure of the binary mixture based on the type 
of percolating cluster, characterised by the largest eigenvalues 
of T, Til and 722. 



III. RESULTS 

The phase diagram of the pure fluids was analysed in 
Ref. [22|. The fluids exhibit a first order phase transi- 
tion that ends at a critical point. The transition involves 
two fluid phases at different densities and fraction of un- 
bonded sites. One phase (high density, small fraction of 
unbonded sites) is always percolated while the other (low 
density and large fraction of unbounded sites) is perco- 
lated only at pressures and densities close to the critical 
point. We call these phases network fluid and vapor, re- 
spectively. 

The results for the binary mixture are presented as 
temperature-composition phase diagrams at constant 
pressure and the following remarks apply to all of them: 
the binodal lines are depicted as solid-black lines; shaded 
areas are two-phase regions; empty circles are critical 
points; black squares are azeotropic points. In all cases 
the network fluid phases are percolated (the liquid side of 
the binodal is percolated). Therefore, at a given compo- 
sition, the mixture is percolated at temperatures below 




Gel Bicontinuous gel Mixed gel 



FIG. 3. Schematic representation of the gel structures: a standard gel of species 1 (left), a bicontinuous gel or bigel (middle), 
and a mixed gel (right). The figure is a cartoon consistent with the bonding rules of the model. Each particle is bonded to three 
other particles (bonding sites not shown) and there are no closed loops. The upper row depicts the position of the particles. 
Bonds between particles of the same species are represented in the lower row by straight lines connecting the center of mass 
of two bonded particles. Each red (blue) line corresponds to a bond between two particles of species 1 (2). Bonds between 
particles of different species are not shown. 



the temperature of the network fluid- vapor phase tran- 
sition. In order to distinguish different percolated struc- 
tures, percolation thresholds for each species are also cal- 
culated. The percolation line of species 1 is depicted as 
a solid-red line. When we remove the bonds between 
species 1 and 2, species 1 is percolated on the right of 
this line (indicated by the red arrow). The dashed-blue 
line is the percolation line of species 2. On the left of this 
line (indicated by the blue arrow), species 2 is percolated 
when we remove the bonds between species 2 and 1. 



A. Symmetric binary mixture AAA-BBB 

We start by describing the results for a simplified 
model: a binary mixture where the particles of species 1 
have three patches of type A and those of species 2 have 
three patches of type B. In addition, we analyse the sym- 
metric case only. That is, we set saa = ^bb = £ and vary 
the bonding energy of the AB bonds {s\^ = sab/^aa)- 

Two representative phase diagrams of this mixture are 
depicted in Fig. |4]for different bonding energies: = 
0.97 (top) and e\j^ = 1.03 (bottom). 

In the first mixture (top panel) the interaction en- 
ergy between sites on different species is lower than the 
interaction energy between sites on the same species 
i^AB — 0-97). At high temperatures there is a network 
fluid-vapor phase transition with a positive azeotrope 
at X = 0.5. At low temperatures a demixing region, 
bounded by an upper critical point, occurs between two 
network fluids. The demixing is driven by the bonding 
energy: at low temperatures most sites are bonded but 



AB bonds hardly occur as they increase the bonding en- 
ergy. Obviously, the same effect occurs at intermediate 
temperatures (below the network fluid- vapor phase tran- 
sition and above the demixing region). However, in this 
range of temperatures the level of association is much 
lower and the energy gain does not compensate the loss 
in the entropy of mixing. As a result, the mixture is 
stable. 

Let us focus now on the percolation threshold(s). Be- 
low the network fluid- vapor phase transition the mixture 
is percolated. The analysis of the percolation threshold 
for each species (solid-red and dashed-blue lines) reveals 
three different gel structures. In the region located to 
the left of the solid-red line {x ^ 0) only the particles of 
species 2 are percolated. The gel is a standard (G2)- The 
same behaviour occurs (recall that the mixture is sym- 
metric) to the right of the shaded-blue line {x ^ 1), but 
in this case the gel is a standard (Gi). In the interme- 
diate region, both species percolate separately, there is 
a bicontinuous gel (BG). As the temperature decreases, 
the BG region spans a wider range of compositions. At 
very low temperatures, however, the bigel is pre-empted 
by the demixing region where two standard gels Gi and 
G2 coexist. 

The second mixture. Fig. |4] (bottom panel), illustrates 
the behaviour of a compound forming mixture: the at- 
traction between sites on different species is stronger than 
that between sites on the same species. This has a pro- 
found effect both on the thermodynamic and the perco- 
lation behaviour. The network fluid- vapor phase transi- 
tion and the azeotropic point at x = 0.5 are still present. 
However, this is a negative azeotrope, characteristic of 
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FIG. 4. Reduced temperature- composition phase diagram at 
constant pressure for a symmetric AAA — BBB binary mix- 
ture with s'^B — 0-97 (top) and s'Xb — 1-03 (bottom). The 
pressure is p* = pVs/sAA = 4.19 x 10~^, well below the criti- 
cal pressure of the pure fluids, x is the composition of species 
1 (HSs with three patches of type A). See the beginning of 
section IIIII for a description of the symbols and the graphical 
codes. 



FIG. 5. Bonding energies between different sites of the model. 



sible for the connectivity of the spanning cluster. As 
the composition of the mixture increases, the fraction of 
bonds between particles of species 1 {AA bonds) also in- 
creases with a corresponding decrease in the fraction of 
AB bonds. As a result the MG structure is replaced by 
a standard Gi gel. The symmetric behaviour (formation 
of a G2 gel) occurs at low values of the composition. 

The percolation thresholds of species 1 and 2 tend 
asymptotically to x = 2/3 and x = 1/3 respectively, as 
the temperature vanishes. For mixtures with composi- 
tions x<l/3orx>2/3 there are no MG structures, 
as the fraction of AB bonds is small and the gel is a 
standard single species structure. This is not the case 
for bicontinuous gels. At low temperatures the percola- 
tion thresholds of species 1 and 2 tend asymptotically to 
X = 1 and X = (note that in Fig. [4l(top) the percolation 
lines are not shown at temperatures below the demixing 
region). If there was no demixing, a bicontinuous gel 
would be found at any composition at low temperatures. 

The symmetric mixture analysed above is the simplest 
system where bicontinuous and mixed gels are formed. 
However, it is not possible to find both structures in 
the same mixture. In addition, the bicontinuous struc- 
ture is pre-empted by a demixing transition between two 
standard gels. To address these questions, we analyse a 
(more) general class of AAC — BBD mixtures in the next 
section. 



compound forming mixtures. Another important dif- 
ference is that the mixture is always stable at temper- 
atures below the network fluid-vapor phase transition. 
The demixing region disappears as there is no driving 
force for phase separation. The bonding free energy (en- 
tropy of mixing) is minimal (maximal) when the number 
of AB bonds is maximal. 

At the percolation level there are also important differ- 
ences. In this mixture, there is no bigel structure. Near 
X = 0.5 neither species percolates separately although 
the mixture is percolated. This is the mixed gel (MG) 
structure where the AB bonds dominate being respon- 



B. Symmetric binary mixture AAC-BBD 

There are up to 10 bonding energies in AAC — BBD 
binary mixtures and, therefore, many possibilities to de- 
scribe the formation of bicontinuous and mixed gels. We 
have focused on the simplest by restricting the bonding 
energies as depicted in Fig. [5l The bonding interac- 
tion between any pair of sites on particles of the same 
species is identical {eaa = ^ac = ^cc = £bb = £bd = 
£dd = £)• This sets the energy scale of the model; 
The strongest interaction is that between sites of type 
C and D {scd = l.OSe); There is no interaction be- 
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FIG. 6. Reduced temperature- composition phase diagram at 
constant pressure for an AAC — BBD binary mixture with 
— 0.83. The pressure is p* = 4.19 X 10"^ weh below the 
critical pressure of the pure fluids, x is the molar fraction for 
species 1. 

tween sites of type A and D nor between sites B and 
C {sad = ^BC = 0); Finally, we vary sab^ the interac- 
tion between sites of type A and B. The physical idea 
behind this choice is the following: if the interaction be- 
tween unlike species is due to the strongest CD bonds 
only {sab = 0), a bicontinuous gel structure is expected 
at low temperatures. The reason is that most of the par- 
ticles will be bonded to two other particles of the same 
species and to one particle of the other species, result- 
ing in long interconnected chains of identical particles. 
On the other hand, if the interaction between A and B 
sites is sufficiently strong a mixed gel is expected as many 
particles are bonded to at least two particles of the other 
species (through CD and AB bonds). 

In what follows we describe temperature-composition 
phase diagrams for different values of the interaction be- 
tween the A and B sites. The pressure is p* = pvs/s = 
4.19 X 10~^, the same as in the previous section (i.e. well 
below the critical pressure of the pure fluids). 

1- Sab < e 

In Fig. [6] we illustrate the phase diagram for a binary 
mixture with = sab/^ = 0.83. The behaviour is 
qualitatively the same for lower values of the AB inter- 
action down to = 0. At high temperatures and large 
(small) values of the composition there is a phase transi- 
tion between a vapor and a gel phase where only species 
1 (2) is percolated. The regions of stability of Gi and 
G2 gels are relatively small due to the presence of a large 
demixing region which meets the fluid-vapor binodal at 



a G1VG2 triple point. As the temperature decreases, 
the two-phase region increases until it spans a very large 
range of compositions. However, at sufficiently low tem- 
peratures, a new region of stability appears. Based on the 
percolation analysis (the percolation lines are not shown 
as they are always inside the demixing region) this new 
low temperature phase is a bicontinuous gel. The phase 
diagram exhibits a first order phase transition between 
Gi or G2 and the BC phase. 

The phase behaviour can be rationalized in terms of the 
competition between the entropy of bonding and the en- 
ergy of bonding. At very low temperatures, the energy of 
bonding dominates. The mixture is stable due to the for- 
mation of CD bonds which are energetically favourable. 
The stable phase is a bicontinuous gel since each particle 
is bonded to a single particle of the other species. The 
other two bonding sites (type A for species 1 or type B 
for species 2) are bonded to identical sites, on particles of 
the same species, giving rise to two interconnected span- 
ning networks. The reason is that A A and BB bonds are 
stronger than AB bonds. Strictly speaking, the previous 
argument is correct only at x = 0.5 and very low tem- 
peratures. If X > 0.5, for example, it is also possible to 
find bonds of type CC or AC since there are more sites 
of type C than sites of type D. However, it is only at 
values of the composition very close to those of the pure 
fluids that the BC phase is replaced by standard Ci or 
G2 gels. At high temperatures, the behaviour of the sys- 
tem is dominated by the entropy of bonding. CD bonds 
are energetically favourable, but the formation of these 
bonds decreases the entropy of bonding since the other 
two bonding sites are forced to form bonds with identi- 
cal sites {AB bonds are still unfavourable). Note that 
the same argument applies at low temperatures, but at 
high temperatures the gain in the energy of bonding from 
CD bonds does not compensate the loss in the entropy 
of bonding (in addition, the fraction of unbonded sites 
is higher at high temperatures). The entropy of bonding 
drives the phase separation between the two standard 
gels Ci and C2 where bonds between sites of identical 
particles are maximized, independently of their type {e.g. 
AA, AC or CC bonds for the Ci phase). 

As expected, the demixing region decreases as the 
interaction between the A and B sites increases. At 
^AB ~ 0-^^ (^^^ -^^S- [3top) the demixing region is still 
connected to the network fluid-vapor binodal. At inter- 
mediate temperatures there is a triple point where three 
bicontinuous gels at different compositions coexist. Be- 
low the triple point there are two symmetric and small re- 
gions of phase coexistence of two BGs, which end at lower 
critical points. At = 0.89 (see Fig. [71 middle) the 
demixing region is a closed loop of immiscibility bounded 
by upper and lower critical points. At =0.93 (see 
Fig. [3 bottom) the demixing region has completely dis- 
appeared. There is no phase separation at temperatures 
below the network fluid-vapor binodal. 

Let us focus now on the percolation threshold(s). In 
all cases depicted in Fig. [3 the percolation lines inter- 
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FIG. 7. Caption as in Fig. [6] for mixtures with: = 0-87 
(top), = 0-89 (middle), and = 0-93 (bottom). 



cept the network fluid-vapor binodal at high tempera- 
tures and tend asymptotically to x = and x = 1 as 
the temperature vanishes. It is possible to flnd a BG 
at any composition if the temperature is sufficiently low. 
At compositions near x = and x = 1 (pure fluids) the 
stable phase is a standard gel, Gi or G2, where only one 
species percolates. These regions grow as the interaction 
between the A and B sites increases. As a result, the 
intermediate region (where the BG is stable) decreases. 
This is to be expected as the fraction of particles with 
more than one bond between different species increases 
as Sab increases. 

Considering only the percolation thresholds, we con- 
clude that the stability of the bicontinuous gel increases 
as the interaction between the A and B sites decreases. 
However, at low eab there is a two standard gels demix- 
ing region, which effectively reduces the stability of the 
BG phase. Consequently, there is an optimal value of 
€ab that maximizes the thermodynamic stability of the 
BG structure. 



2. Sab ~ s 

The most interesting systems correspond to mixtures 
where the energy of the AB bonds (between unlike 
species) is similar to the energy of the AA and BB bonds 
(between like species). In Fig. [8] we illustrate the phase 
diagram of a binary mixture with = 0.99 (top) and 
= 1.02 (bottom). In both cases we flnd the four dif- 
ferent gel structures by varying the temperature and/or 
the composition of the mixture. In the pure fluid regime 
(x ^ and x ^ 1) only one species is percolated giv- 
ing rise to standard Gi or G2 gels, while in the equimo- 
lar regime around x = 0.5 there is a competition be- 
tween bicontinuous and mixed gels. The bigel is stable 
at high temperatures and the mixed gel at low tempera- 
tures. The stability of the BG or MG is again the result 
of a competition between the entropy and the energy of 
bonding. 

Let us start by describing the mixture with = 1.02, 
that is, all bonds between different species {AB and CD) 
are stronger than the bonds between particles of the 
same species {e.g. A A, AG or GG). We expect a stable 
mixed gel (or standard gels if the composition is far from 
X = 0.5). In fact, a MG appears at low temperatures, 
but at high temperatures there is a small region where a 
BG is stable. The stability of the BG results from a gain 
in the entropy of bonding. Consider a particle of species 
1. It has 2 sites of type A which form AA, AG and AB 
bonds, and 1 site of type G which forms CC, GA or CD 
bonds. There are 27 different combinations when all sites 
are bonded. 7 of these favor the formation of mixed gels 
(those with 2 or 3 bonds between different species) while 
20 favour bicontinuous gels (those with no more than 1 
bond between different species). At high temperatures, 
the entropic part of the bonding free energy dominates 
and it is possible to stabilize bicontinuous gels even when 
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FIG. 8. Caption as in Fig. [6] for mixtures with: e*AB — 0-99 
(top), Sab = 1-02 (bottom). 



the AB bonds are stronger than, for example, AA or AC 
bonds. The opposite behaviour (a stable mixed gel when 
the AB bonds are weaker than the AA or AC bonds) is 
also possible. In Fig. [8] (top), we illustrate the phase 
diagram of a mixture with =0.99. At high temper- 
atures the BG is stable but at low temperatures we find 
a small region of stability of the MG. Finally, at zero 
temperature, the BG reappears at x = 0.5 (note that 
the ground state, at x = 0.5, is a bicontinuous gel since 
there are only CD and AA or BB bonds in order to min- 
imize the energy). The MG is stable in a small region 
at low temperatures because it minimizes the energy of 
bonding. Recall that there are 7 combinations of bonds 
that favour a MG. When e 



AB 



0.99 the average energy 
per bond of these configurations is {s\^q) ~ 1.007, while 
the average over the configurations that favor a BG is 
{ela) = 1.002. 
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FIG. 9. Reduced temperature-composition phase diagram at 
constant pressure for an AAC — BBD binary mixture with 
^AB = 1-02. The pressure is p* —pvs/e = 5.24x10"^, slightly 
below the critical pressure of the pure fluids. 



If the above argument is correct, we should find a sta- 
ble MG at low temperatures when = 0.975 (at this 
value of Sab the mixed and bigel average energies per 
bond are equal, {s^^q) = {smq) ~ l)- However, the 
MG appears when ^ 0.977 (at the given pressure). 
This discrepancy may be traced to the assumption that 
all the sites are bonded, which is only an approximation 
at finite T. In fact, the fraction of unbonded sites can 
play a major role in the stability of percolated structures. 
For example, at high temperatures and pressures slightly 
below the critical pressure of the pure fluids, the frac- 
tion of unbonded sites is relatively high. Then, a signif- 
icant fraction of particles will have only one or two sites 
bonded, increasing the number of configurations that fa- 
vor a mixed gel. As a result the MG becomes stable 
at temperatures above the region of stability of the BG. 
An example is shown in Fig. [9l where the interaction 
between sites A and B is = 1.02 (the same as in 
Fig. [8] (bottom)) at a pressure p* = pVs/e = 5.24 x 10~^ 
(slightly below the critical pressure of the pure fluids). 
Lowering the temperature at x = 0.5 we find a reentrant 
sequence of gel structures: MG — BG — MG. 

In Fig [To] we represent a different phase diagram for the 
mixture with = 1.02. This is a cut of the (T^rj^x) 
phase diagram at constant composition x = 0.47 (note 
that in this representation a tie line connecting two co- 
existing points is out of the plane). The dotted-black line 
is the percolation threshold of the mixture (not shown in 
the other phase diagrams since the fluid phase is always 
percolated). 
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FIG. 10. Reduced temperature-packing fraction phase dia- 
gram at constant pressure of the symmetric AAC — BBD 
binary mixture at composition x = 0.47. The interaction be- 
tween the A and B sites is s^b — 1-02. A tie hne connecting 
two coexisting phases is out of this plane. 
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FIG. 11. Caption as in Fig. [6] for mixtures with: s^b = 1-05. 



3. £AB > £ 

When the interaction between sites A and B is much 
stronger than the interaction between sites of the same 
species, we expect only mixed and standard gels. We fin- 
ish this section with a simple case where = = 
1.05, that is, all the bonds between different species 
are stronger than the bonds between identical particles. 
This system is similar to the symmetric binary mixture 
AAA — BBB^ analysed in section IIII A| but not identi- 
cal. The phase diagram is depicted in Fig. [TTJ There 
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X 

FIG. 12. Reduced temperature-composition phase diagram 
at constant pressure (p* = 4.19 x 10~^) of an asymmetric 
AAC — BBD binary mixture. See the text for details of the 
bonding interactions. 



are 3 azeotropic points on the network fluid-vapor bin- 
odal (2 positive azeotropes at x ^ 0.18 and x ^ 0.82, 
and a negative azeotrope at x = 0.5). The mixture is 
completely stable below the network fluid- vapor binodal 
temperature. A mixed gel is stable in a broad region 
around x = 0.5. This region is bounded on the left by 
a G2 structure and on the right by a Gi structure. As 
expected, no BG structure is found. 



C. Asymmetric binary mixture AAC-BBD 

The phenomenology described above is not restricted 
to symmetric binary mixtures, which were considered for 
simplicity only. In Fig. [12] we show an example of a 
similar phase diagram for an asymmetric binary mix- 
ture. It corresponds to a mixture where the symmetry 
is broken through the interactions between the sites on 
particles of species 2. The set of interaction energies is: 
£aa = £ac = ^cc = £ (it sets the scale of energy), 
£bb = £bd = £dd = 0.75£, scd = 0.95^, eab = 0.886 
and Sad = £bc = 0. The mixture is always percolated 
at temperatures below the network fluid-vapor binodal 
(which now occupies a significant region of the phase di- 
agram). By varying the temperature and composition 
of the mixture it is possible to find the four percolated 
structures described for symmetric mixtures. As before, 
the stability of the different structures can be understood 
in terms of the competition between the energy and the 
entropy of bonding. 
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IV. CONCLUSIONS 



We have studied the thermodynamics and percolation 
thresholds of a simple model of patchy colloidal particles: 
a binary mixture of hard spheres with different types of 
bonding sites on their surfaces. Despite the simplicity of 
the model, the mixture exhibits a very rich phase dia- 
gram. The network fluid regime includes up to four gel 
phases, including bicontinuous and mixed gels. The bi- 
continuous gels described here are similar to those stud- 
ied by Hall an co-workers (26l . [27l | in mixtures of dipolar 
colloids. 

The stability of the different percolated structures is 
determined by a competition between the entropy of 
bonding (number and type of bonding sites) and the en- 
ergy of bonding (bonding energy). Therefore, we expect 
the results to be relevant to a wide range of patchy par- 
ticle systems, regardless of the size and/or the geometry 
of the colloids. 

We have shown that by tuning the bonding energies it 
is possible to design binary mixtures with stable bigels 
and/or mixed gels in a wide range of temperatures and 
compositions. When the strength of the bonds between 
like and unlike particles is similar, we find an interesting 
competition between the mixed and bicontinuous gels, 
including reentrant behaviour. The mixtures where both 



bicontinuous and mixed gels compete are promising can- 
didates to fabricate materials with novel physical prop- 
erties. An obvious example is a mixture where one of 
the species transmits a given property (light, electric- 
ity...) while the other does not. It is possible to control 
the transmission of that property by simply varying the 
temperature, in the regime where MG and BG compete. 

We have not considered here the stability of position- 
ally ordered phases. Solid phases will appear at low tem- 
peratures and/or high pressures preempting part of the 
phase diagrams. Nevertheless, the gels we have analysed 
appear at relatively low packing fractions. In addition, 
it is expected that the random distribution of patches on 
the particle surfaces will frustrate the formation of solid 
structures. 
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